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Abstract 

Five years ago, Lane Hughston showed that some of the symmetric 
informationally complete positive operator valued measures (SICs) in di- 
mension 3 coincide with the Hesse configuration (a structure well known 
to algebraic geometers, which arises from the torsion points of a certain 
elliptic curve). This connection with elliptic curves is signalled by the 
presence of linear dependencies among the SIC vectors. Here we look for 
analogous connections between SICs and algebraic geometry by perform- 
ing computer searches for linear dependencies in higher dimensional SICs. 
We prove that linear dependencies will always emerge in Weyl-Heisenberg 
orbits when the fiducial vector lies in a certain subspace of an order 3 
unitary matrix. This includes SICs when the dimension is divisible by 3 
or equal to 8 mod 9. We examine the linear dependencies in dimension 
6 in detail and show that smaller dimensional SICs are embedded in this 
structure, potentially impacting the SIC existence problem. We extend 
our results to unitary matrices that are not order 3. Finally, we align our 
work with recent studies on representations of the Clifford group. 
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1 Introduction 



Symmetric informationally-complete positive operator-valued measures (SICs) 
[TJ[2] represent a general form of measurement in quantum theory. As the 
more familiar projective measurements are associated to an orthonormal ba- 
sis in Hilbert space, a SIC is associated to an over-complete set of N 2 unit 
vectors in T-L N such that the scalar product between any distinct two is always 
constant, i.e. 



SICs have practical applications in quantum state tomography [3,4 , quantum 
communication 0H], quantum cryptography classical high precision radar 
[9, 10J and classical speech recognition [TT]. A significant amount of work is 
aimed towards proving their existence in all dimensions although no general 
proof is currently known. Numerical studies have successfully found SICs when 
N < 67 [12] and analytical solutions are known for 20 of these dimensions (see 
references in [12] plus one recent solution for N — 16 in [13J). Such results 
promote the belief that SICs can always be found, but the current solutions are 
rather dimension-dependent and do not provide an overall coherent picture. 

All known SICs are group covariant, meaning they can be obtained from 
the action of a group on a single fiducial vector. The vast majority of SICs are 
covariant with respect to the Weyl-Heisenberg (WH) group and if the dimension 
is a prime, it has been shown this is the only possible group [TJ] . 

There are a few curious properties of SICs that are not yet understood, but 
clearly hint at an underlying mathematical structure. This scenario could sup- 
port suggestions that SICs should play a deeper role in quantum mechanics [15] . 
The most well-known of these properties is the conjecture of Zauner symmetry, 
which states that every SIC vector is an eigenvector of a certain order three uni- 
tary matrix jTl fTSHTB] . Another property stems from current SIC solutions; in a 
preferred basis, the components of the SIC vectors are expressed as nested rad- 
icals. Such complex numbers are quite exceptional and arise because the Galois 
group of the polynomial equation one must solve is solvable |12U17| . However, it 
is a third observation that we shall focus on in this paper, namely the connection 
between SICs and algebraic curves in dimension 3. 

The WH group appears naturally in the theory of elliptic curves [18 , sug- 
gesting this might be an advantageous framework for the study of SICs. Further- 
more, there is a very special connection between elliptic curves and linear depen- 
dencies among SIC vectors in dimension 3, first noted by Lane Hughston |19| . 
This leads to the Hesse configuration in the complex projective plane — a partic- 
ular structure well-known in algebraic geometry. We aim to explore this special 
connection, looking for analogous linear dependencies in higher dimensions. It 
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was our hope that another configuration in complex projective space might re- 
veal some new structure behind SICs and broaden the areas from which they 
can be studied. This did not turn out quite as planned for all the SICs we looked 
at but we do see similarities for certain SICs in dimension 8, which could sig- 
nal another connection to elliptic curves. No sets of N linearly dependent SIC 
vectors exist for N = 4. 5 or 7, but we do find them again when N — 6, 9 and 
12 — however, in these latter cases the linear dependencies are not uniquely tied 
to SICs. Additionally, we found an unexpected method of generating 2- and 
3-dimensional SICs from the linear dependencies among vectors in 6- and 9- 
dimensional SICs, respectively. We don't know whether this can be generalised 
to all dimensions divisible by 3. 

In section 2 we review some basic facts about the Weyl-Heisenberg and 
Clifford groups, and explain notations and terminology used in the rest of the 
paper. Section 3 introduces the Hesse configuration and linear dependencies 
in SICs in dimension 3. Section 4 proves that linear dependencies can always 
be found in WH orbits when the fiducial vector is an eigenvector of an order 3 
unitary matrix. In dimensions divisible by 3 and dimensions equal to 8 mod 9, 
these WH orbits include SICs. In section 5, we give numerical results from our 
computer search for dependencies in dimension 4 to 8. We focus specifically on 
data relating to the SIC case in dimensions 6 and 9. In section 6, we prove the 
sets of linearly dependent vectors in dimension 6 contain 2-dimensional SICs 
and show this also occurs in dimension 9. We do not have data for dimension 
12 or higher, but see no reason why this pattern shouldn't continue for all 
dimensions divisible by 3. In which case, this could open a new avenue for a 
SIC extistence proof. Section 7 turns away from SICs altogether and proves 
that linear dependencies will always arise in WH orbits with a fiducial vector 
invariant under a unitary matrix whose order is a divisor of N in dimensions 
where N is a product of distinct prime numbers. We give an example of this 
using linear dependencies in even dimensions. Section 8 relates our results to 
a recent representation of the Clifford group |13| . Section 9 summarises our 
results and points at some possible clues for the SIC existence problem. 



2 The Weyl-Heisenberg and Clifford groups 

In this section we review some basic facts about the WH and Clifford groups 
which will be needed in the sequel. For more details see (for example) [16] . Let 
|0), . . . \N — 1) be the standard basis in dimension N. Define the operators X 
and Z by 

X|«) = |« + l> (3) 
Z\u) = uj u \u) (4) 

where addition of ket-labels is mod N and u> = e^. The Weyl-Heisenberg 
displacement operators are then defined by 

D p = T P1P2 X P1 Z P2 (5) 
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where p = (pi), t = — e « . With this definition 

D p £» q =r^Z> p+q D l = D - P (6) 

where the symplectic form (•, •) is defined by 

(P,q) =P2qi -Pi«2 (7) 
It is convenient to introduce the notation 

N=i N N ° dd (8) 

[27V 27V even v ' 

We then have t n = 1 and 

D p = £> q if p = q mod N (9) 

(note that if N is even t n = — 1 and p = q mod N implies D p = ±_D q ). We 
define the Weyl-Heisenberg group to consist of all operators of the form T n D p , 
with n E Zfj, p G Z 2 ^ (where Z^y is the set {0, 1, . . . ,N — 1} equipped with 
arithmetic modulo N). 

In many situations the phase factor in the first of Equations © is not im- 
portant. It is accordingly convenient to introduce the notation A=B to signify 
that A and B are equal up to a phase. We then have 

DpD^D^Dp^Dp+q (10) 

for all p. q. Neglecting phases in this way the WH group reduces to the iV 2 op- 
erators {Dp : p g Z 2 ^} (in more technical language these facts can be expressed 
by the statement that the collineation group is = Zjv x Zjv). Similarly, if we 
ignore phases then the orbit of a vector \ip) under the WH group reduces to the 
set {D p \ijj) : p e Z 2 ^} (in more technical language this set is the orbit of \tp) 
considered as an element of (DP^ -1 under the action of the collineation group). 
In the sequel, by an abuse of terminology, we shall refer to this set as the WH 
orbit of \ip). 

The symplectic group SL(2, Z^) consists of all matrices 

c =(°*) (ii > 

such that a, j3, 7, S G Zjy and detG = 1 (mod N). To each such matrix there 
corresponds a unitary Ug, unique up to a phase, such that 

UgDpUl = D gp (12) 
If p is relatively prime to TV we have the explicit formula 

Ug = ^Yl r t3 "( Su2 - 2uv+av2 )\u){v\ (13) 

* u,v— 
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where /3 _1 is the multiplicative inverse of (3 (mod N) and e t9 is an arbitrary 
phase. If (3 is not relatively prime to N we use the fact [16 that Q has the 
decomposition 

i 2 ) ("> 

where /?2 both are relatively prime to N so that {7g 15 Ug 2 can be calculated 
using Equation (fT3j) . J7g is then given by 

Ug = U gi Ug 2 (15) 

up to an arbitrary phase. We shall refer to unitaries of the form e* Ug as 
symplectic unitaries. The Clifford group then consists of all products e t6 D p Ug. 
We have 

(D p Ug) (D p ,Ug,)=D p+gp ,Ugg, {D p Ug) ] =D_g-^Ug-, (16) 

for all p, p\ 

As mentioned in the Introduction it turns out that every known WH SIC 
fiducial vector is an eigenvector of a canonical order 3 Clifford unitary and, 
conversely, that in every dimension where an exhaustive search has been made 
every canonical order 3 Clifford unitary has a SIC fiducial vector as one of its 
eigenvectors |12II16| . In these statements the term "canonical order 3 Clifford 
unitary" refers to a unitary of the form 

U = e lB D p Ug (17) 

where p is arbitrary and Q is any element of SL(2,Zjy) with the property 

Tr(£) = -1 mod N (18) 

If N = 3 one needs to impose the additional requirement that Q ^ I. It is 
then guaranteed that, with a suitable choice of the phase e l9 , U is order 3. The 
matrix Q is always conjugate |12II20| . either to the Zauner matrix 

«-(? :0 w> 

or, in the case of dimensions N — 9k + 3, with k > 1, to the matrix 

A ={~N N N k) ( 2 °) 

In this paper we will confine ourselves to the case of fiducials which are eigen- 
vectors of Uz- Choosing the phase in Equation (113]) to be e = e 12 one 
finds [T] that the eigenvalues of Uz are l,rj,rj 2 , where rj = . We denote the 
three corresponding eigenspaces as "Hi,?^, and H^, and give their dimensions 
in Table 1. Except in dimensions equal to 8 mod 9, fiducials are only found 
in the highest dimensional eigenspace(s) of Uz- When N = 8 mod 9 fiducials 
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Table 1: Multiplicities of the eigenvalues of Uz for different dimensions. 



are found in all 3 eigenspaces (T^j. We will hereafter refer to (which, the 
numerical data suggests, always contains SIC fiducials) as the Zauner subspace. 

Finally, let us note that Z has non-trivial fixed points if and only if the 
dimension is divisible by 3. Specifically: 



Zp 




N ^ mod 3 



N = mod 3 



(21) 



3 Linear dependencies in dimension 3 



Lane Hughston showed that vectors in some 3-dimcnsional SICs coincide ex- 
actly with the nine inflection points of a particular family of elliptic curves |19| . 
The pattern of linear dependencies among the SIC vectors reproduces the Hesse 
configuration [21 , 22J and so we shall first look in more detail at linear depen- 
dencies coming from SICs in dimension 3. This section forms our motivation 
for studying higher dimensional dependency structures and, although SICs in 
dimension 3 have been studied fairly extensively before now |14[[T6]. we hope 
the following approach provides a useful way of looking at things. 

In dimension 3, there is a continuous one-parameter family of SICs. The 
nine explicit (un-normalised) vectors are 
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The parameter 6 needs only be considered in the interval [0, ^] as all other 
SICs are equivalent to these via transformations in the extended Clifford group. 
It is clear from inspection that any SIC in dimension 3 will have three sets of 
three linearly dependent vectors. However, there are additional dependencies 
for certain choices of 9. We can find these values by setting the determinant of 
any three other SIC vectors to zero. This gives 
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and we find the condition 
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€{0,1,2} (24) 




Given our limits for 9, we are left with the options 8 — or in order to obtain 
a 'special' SIC with 12 sets of three linearly dependent vectors. The SIC with 
9 = is a particularly special case as its fiducial vector 



(25) 



is invariant under the entire symplectic group. We can see this by rewriting the 
projector onto this vector as 

\<f>}(<f>\ = l-U v V= ( - 1 \ ) (26) 

The matrix V is called the parity matrix, and its unitary representative U-p is one 
of the phase-point operators introduced by Wootters (23] . It is unique in being 
the only element (other than the identity element) invariant under conjugation 
by every other element in the symplectic group and so forms the group's centre. 
As the right hand side of Equation (|26p is invariant under SL(2, Zjv), the SIC 
fiducial vector \(j)) must be an eigenvector of every symplectic unitary matrix. 
Of interest to us is its invariance under the action of the symplectic canonical 
order 3 unitaries. There are eight such unitaries altogether (including the Zauner 
unitary) but we can consider the four unitaries 
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since the other four are squares of these. Note that Ui o 2 \ corresponds to 

\1 2 ) 

the Zauner matrix in Equation (fT9|) of Section 2. The SIC fiducial defined 
in Equation is then precisely the intersection point of the 2-dimensional 
subspaces Hi from each of these four canonical order 3 unitaries. We can 
then ask whether any other SIC vectors are also contained within these same 
subspaces. It turns out that the remaining SIC vectors are distributed evenly 
amongst the subspaces: 2 further vectors per unitary subspace. 

When we view the full SIC 3 lattice in the complex projective plane 

(where each point is a SIC vector and each line is a 2-dimensional subspace) we 
generate the Hesse configuration. In Figure [T] the fiducial vector is found at 
the point at the bottom left of the plane. The left-hand image shows the three 
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linear dependencies in every dimension 3 SIC while the right-hand image shows 
the pattern of symplectic canonical order 3 unitary subspaces that intersect 
at the fiducial vector \<fr), i.e. the SIC fiducial for which 9 — 0. The Clifford 
group also contains WH conjugates of the form D p UD^ p , where U is one of the 
four unitaries in Equation (J28[) , so in this way we find a total of 12 order three 
unitaries and hence 12 subspaces. Remarkably, they always intersect in sets of 
four at a point described precisely by one of the SIC vectors. We recover the 
famous Hesse configuration of 9 points and 12 lines in the projective plane |21| . 
To obtain this full picture from Figure Q] we must translate the lines in the 
right-hand image. Each line can be translated three times, so we find 12 lines 
in total. 




Figure 1: A SIC in the complex projective plane. The left-hand image shows the 
three sets of three linear dependencies in every SIC and the right-hand image 
shows the four subspaces from Equation (j28|) that intersect at the fiducial vector 
with 9 = 

Hesse originally discovered his configuration by examining elliptic curves. In 
particular, the family of cubics 

x 3 + y 3 + z 3 + Xxyz = A e <D (28) 

is invariant under the WH group and the determinant of their Hessian (the 
matrix containing the second derivatives of the curve) is also a cubic in the 
same family. Bezout's theorem states that two cubics in (DP 2 will intersect 
in nine points and these nine inflection points of the cubic in Equation (|28p 
coincide exactly with the nine SIC vectors. There are four singular members in 
this family, known as the Hesse pencil, that arise when A = oo and A 3 = — 1. 
Each singular curve, Ti, degenerates into three projective lines 

To : xyz = 

Ti : (x + y + z)(x + r\y + t] 2 z){x + rfy + rjz) = , 2Q . 
T2 : (x + rjy + z)(x + rj 2 y + rj 2 z){x + y + rjz) = 
T 3 : {x + i] 2 y + z)(x + rjy + ijz){x + y + rfz) = 

Each line passes through three inflection points and each inflection point lies on 
four lines, forming our familiar Hesse configuration. It is often denoted (94, 123) 
to reflect the property that the 9 points each lie in 4 lines and the 12 lines each 
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contain 3 points. A more detailed survey of the Hesse configuration is given 
in [22] and its connection to SICs in [24] . 

Returning to the symplectic canonical order 3 unitaries, this pattern of sub- 
space intersections has a connection to mutually unbiased bases. Each of the 
twelve lines in the Hesse configuration uniquely determines an orthogonal vector 
up to a scalar, which we shall call a normal vector from here onwards. These 12 
normal vectors coincide with the standard MUB in dimension 3, obtained from 
the eigenvectors of the cyclic subgroups of the WH group. Note that we use 
the acronym MUB to mean a complete set of N + 1 mutually unbiased bases. 
As the MUB vectors lie orthogonally to the 2-dimensional symplectic canonical 
order 3 unitary subspaces, they must sit in the 1-dimensional subspace W^. 

So far we have been looking at the single SIC with 9 = 0. However, one 
also gets a Hesse configuration for each of the 8 SICs with 9 = . Most SICs 
in dimension 3 are part of an extended Clifford group orbit containing eight 
SICs |I6]. The exceptions are the SICs with 9 = and 9 = ?f in Equation (|22|) . 
which lie on orbits containing one and four SICs, respectively. The 'special 
SIC, with 9 = is then on an extended Clifford orbit with seven other SICs 
and, in fact, all eight of these SICs are special in the sense of having 12 linearly 
dependent sets of three vectors. This extends to the MUB connection. Above, 
we found a MUB from the normal vectors — lines in the Hesse configuration — of 
linear dependencies in the 'special SIC with 9 = 0. Note that this MUB is itself 
a single orbit of the extended Clifford group. We can extract eight more MUBs 
in the same way from the other eight 'special SICs'. They can be arrived at 
in a different way. The Clifford group contains altogether 36 order 3 subgroup 
with elements of the form D p U, with U a symplectic order 3 unitary. So far 
we have used only 12 of these. The remaining 24 have non-degenerate spectra, 
and their eigenbases can be collected into 8 sets of 3 mutually unbiased bases. 
Each such triplet is in itself an orbit under the WH group, and forms a MUB 
when a basis from the standard MUB is added |25|[26] , This ties in exactly 
with the eight 'special SICs' with 9 = ^ through the Hesse configuration [25] . 
Each SIC, together with one of these WH orbit MUBs, reproduces the pattern 
(94, 123). It is perhaps interesting to note that we could have arrived at the 
nine 'special SICs' from knowledge of the symplectic canonical order 3 unitaries 
and the additional WH orbit MUBs alone. 

4 Linear dependencies from the eigenvectors of 

U z 

In this section we show how linear dependencies will arise in every dimension. 
The results in this section apply to all orbits in which the fiducial vector lies 
in certain eigenspaces of Uz, depending on the dimension, and not only to SIC 
fiducials. 

From section 2, we can label a vector in a WH orbit, up to a phase, by a 
2-component vector p in 7? N . In this representation, we are interested in the 
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orbit of p under the action of Z. If p is one of the fixed points of Z it will 
form a singlet, otherwise {p, Zp,Z 2 p} forms a triplet since Z is of order 3. If 
the dimension N is divisible by 3, Z has exactly 3 singlets, which are given by 
Equation (|2ip . In other dimensions, the only singlet is p = (0, 0). In the Hilbert 
space, we will use the same terminology to call D p a singlet if p is a singlet, 
and to call {D p \ip) , Dz P \ip) , Dz 2 p |^)} a triplet otherwise. 

Theorem 1. In dimension N = 3fc, any subset of N vectors in a WH orbit 
whose fiducial vector is an eigenvector of Uz is linearly dependent if it contains 
k triplets, or k-1 triplets and 3 singlets. 

Proof. Let N = 3k, and let Uz \ipo) = A IV'o)- When p is a singlet, D p \ipo) lies 
in the same eigenspace of Uz as j^o) because 

U Z D P |^o) = U z D p ulUz \iM = \D Zp |Vo) = A£ p |Vo) (30) 

When p is in a triplet, we construct three new linear combinations of vectors in 
the triplet {D p , Dz P \ip) , Dz 2 p Our vectors are 

\r) = ,D p |V;o) + U z D p |^o) + U 2 z D p |Vo> 

|s) = D p |^o) + V 2 UzD p |Vo) + VU Z D P l^o) (31) 

\t) = D p |^o) + VU Z D P |^ ) + V 2 U 2 Z D P \f ) 

which, it is straightforward to check, are evenly distributed among the three 
eigenspaces Hi, H v and H^, shown in Table 1. It is also clear that the linear 
span of the vectors {D p \ ipo) , Dz P \ ipo) , D z z p l^o)} equals that of the vectors 
i\r),\s),\t)}. 

If a set of N vectors contains k triplets, this gives fc-many of each r-,s-, and 
t-type vector. From Table 1 we know that the r-type vectors lie in an eigenspace 
of dimension k + 1, the s-type vectors lie in an eigenspace of dimension k, and 
the t- type vectors lie in an eigenspace of dimension k — 1 . It is clear then that 
the k r-type vectors cannot fully span their subspace while the k i-type vectors 
are overcomplete and therefore linearly dependent. 

If a set contains k — 1 triplets and 3 singlets, this gives (k — l)-many of each 
r-,s-, and i-type vector, plus the 3 singlets that are in the same eigenspace as 
the fiducial vector \ipo). Therefore there will be k + 2 vectors lying in the same 
eigenspace. Since the largest eigenspace of Uz has dimensionality k + 1, this 
results in linear dependency and concludes the proof. □ 

When the fiducial vector lies in H v or H 2 we can obtain dependencies using 
fewer than N vectors. For example, if \ipo) € T-L^ then linear dependency will 
also occur when the set contains k — 1 triplets and 2 singlets, or k — 2 triplets 
and 3 singlets. And if l^o) S Hri 2 i linear dependency will also occur when the 
set contains k — 1 triplets and 1 singlet, k — 2 triplets and 2 singlets, or k — 3 
triplets and 3 singlets. 
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Theorem 2. In dimension N — 3k + l, any subset of N vectors in a WH orbit 
whose fiducial vector lies in the eigenspace H v or H v 2 is linearly dependent if it 
contains k triplets and 1 singlet. 

Proof. In dimension N = 3k + 1, the three eigenspaces Hi, H n and Hrf have 
dimensionality k + 1, k, and k, respectively. If the fiducial lies in H v , which is 
the only singlet in this case, then together with k s-type vectors it will give k + 1 
vectors in this subspace, resulting in linear dependency. The same argument 
applies when the fiducial vector lies in H 7 f . □ 

Theorem 3. In dimension N = 3fc + 2, any subset of N — 1 vectors in a WH 
orbit whose fiducial vector lies in the eigenspace H„2 is linearly dependent if it 
contains k triplets and 1 singlet. 

Proof. In dimension N = 3fc + 2, the eigenspace H„z has dimensionality k. The 
k f-type vectors, together with the fiducial vector, form a set of k + 1 vectors in 
H v 2, thus they are linearly dependent. □ 

Lastly, we want to note that if a set is linearly dependent, all other sets in 
its WH orbit are also linearly dependent, because they can be translated to the 
first set by the action of the displacement operators D p 's. In Theorem Q] sets 
that contain fc triplets, or k — 1 triplets and 3 singlets, have the property that 
they remain invariant (as a set) under the action of Uz- We shall refer to these 
sets as Zauner invariant sets, and their orbits as Zauner invariant orbits. 



5 Numerical linear dependencies 

In the previous section, we proved that we can find N linearly dependent vectors 
in a WH orbit when the fiducial vector is in a certain Zauner subspace. However, 
this doesn't quite capture the whole picture. In this section, we provide the 
outcome of our numerical search for linear dependencies in dimensions 4 to 8 
(with partial results when N = 9 and 12) and give additional linear dependencies 
that are not covered in the previous section. We are especially interested in 
dimensions divisible by 3, where dependencies are always found among SIC 
vectors, and we discuss details of linear dependencies in dimensions 6 and 9. 
The additional dependencies in these dimensions do not depend on whether 
the WH orbit is a SIC or not, although we find some interesting orthogonality 
relationships in the dependency structure that are unique to SICs (and hence 
do not occur in WH orbits starting with a non-SIC fidicuial vector in Hi). 
Our investigation into dependencies in dimensions 6 and 9 also led to smaller- 
dimensional SICs, but we will postpone a discussion of these until the next 
section. In dimension 8, there are SIC fiducials in the Zauner subspace H,f and 
our numerical data show that these SICs have more linearly dependent sets of 
8 vectors than WH orbits with a vector in H„2 . Recall that this was precisely 
the situation with the 'special SICs' in dimension 3 and so could indicate a 
dimension 8 analogue of the Hesse configuration. 
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In each dimension N, our computer program started with a vector from one 
of the Zauner subspaces, generated the full orbit under the action of the WH 
group and then performed an exhaustive search for all sets of N vectors that are 
linearly dependent. For each dimension, we used vectors randomly chosen from 
the subspaces of the Zauner unitary. In the case of SIC vectors, we used the 
fiducials given in ; where there was a choice of Clifford orbits, we repeated 
the calculation with fiducials from each orbit. 

The numerical results often showed a higher number of linear dependencies 
than accounted for by Theorems 1-3. Table 2 shows the total number of sets 
of N linearly dependent vectors in WH orbits with fiducial vectors in different 
Zauner subspaces. The number predicted by the proofs in the previous section, 
or an upper bound on this number, is given in brackets. 
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24756984 


(68) 


(4200) 




(5796) 


(<766080) 



Table 2: Linear dependencies in a WH orbit when the initial vector is taken 
from each Zauner subspace. The numbers in brackets are the total number of 
sets (or an upper bound on this number) predicted in section 4. 

Our results show that SICs in dimensions divisible by 3 contain linearly 
dependent vectors. As we are primarily interested in SICs, we turn to look 
in more detail at these dimensions. When N = 6 there are 984 sets of six 
linearly dependent SIC vectors, to a numerical precision of 10~ 15 . All the results 
reported in this section are numerical and hold to within this precision. This 
is a higher number of dependencies than we proved must exist in the previous 
section, which finds only 768 of the 984 sets we generated numerically. Recall 
that the proof required the N vectors in the WH orbit to be invariant under the 
Zauner unitary or one of its WH conjugates; the additional 216 sets in dimension 
6 are not of this form. Curiously, they are instead invariant under an order 6 
unitary matrix Um whose symplectic representation is 

(32) 

Each of the 36 SIC vectors lies in 164 different sets of linearly dependent 
SIC vectors, and clearly each of the 984 linearly dependent sets contains 6 
vectors. In the language of complex projective space there are 36 points and 984 
hyperplanes, forming the balanced configuration (36i64, 984@). If this constitutes 
a known pattern in (DP 5 , we do not recognise it. 
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The 984 linearly dependent sets can themselves be collected into orbits un- 
der the WH group. There are 27 orbits of length 36 and 1 orbit of length 12. 
The short orbit arises because it contains only linearly dependent sets invariant 
under the subgroup {1,X 2 Z 4 ,X 4 Z 2 }. This subgroup commutes with the Za- 
uner unitary defined in Equations (fTSJ) and (|19p . which leaves the fiducial SIC 
vector invariant. Twenty two of these WH orbits contain sets that are invariant 
under the Zauner matrix (or sets invariant under a WH conjug ate D p UzDl of 
the Zauner matrix) ; the other 6 are of course invariant under the action of Um . 

Inspired by the 3-dimensional case, we look at the 984 vectors normal to the 
linearly dependent sets, and perform an exhaustive search for orthogonalities 
between them. There is no basis among them — nor seven MUBs! — but there 
are nine sets of four mutually orthogonal normal vectors. These normal vectors 
all belong to the same orbit under the WH group (of length 36). Each set of 
four mutually orthogonal vectors forms an orbit under the subgroup of the WH 
group that contains order 2 elements only. In addition there are over 20,000 
orthogonal triples of normal vectors, only 712 of which lie in a single WH orbit. 

The resulting structure in dimension 6 is summarised in Table 3, where we 
have labelled each WH orbit from 1 to 28. The first orbit is the short one and 
we see that every set of six linearly dependent vectors in this orbit is formed 
from two orbits under the subgroup {1, X 2 Z 4 , X 4 Z 2 }. The orbit containing the 
mutually orthogonal normal vectors is number 11, and orbits 23 to 28 are the 
ones invariant under Um — they contain sets of six vectors forming just one orbit 
under Uz, rather than two. 



WH orbit 


No. orbits under 


No. orbits under 


No. ON 




Uz 


{1,X 2 Z 4 ,X 4 Z 2 } 


quadruples 


1 


2 


2 





2-10 


2 


1 





11 


2 





9 


12-13 


2 








14-22 


2 








23-28 


1 


1 






Table 3: Properties of WH orbits starting with a vector in T~i\ in dimension 6. 

The six-dimensional case shows some similarities to the three-dimensional 
one, but the pattern is not as striking. We now ask to what extent the results 
depend on the fact that we start from 36 vectors forming a SIC. From the previ- 
ous section, we expect at least 768 of the 984 sets of linearly dependent vectors 
to remain unchanged. To check this, we repeated our calculations starting from 
a fiducial vector chosen at random from the Zauner subspace (i.e. in T-Li). This 
will produce 36 vectors typically not forming a SIC. In the cases we looked at, 
we found 984 sets of six linearly dependent vectors in exactly the same pattern 
as when we started with a SIC, i.e. the same p vectors contributed to the same 
sets. 
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However, there are some subtle differences between SICs and non-SlCs in 
dimension 6, despite the dependencies being the same when the fiducial is taken 
randomly from Hi. The four orthogonal quadruples among the normal vectors 
are present only in the SIC case. An additional 216 orthogonal triples vanish 
when we pick the fiducial at random from Hi] the triples all have one normal 
vector in WH orbit 11, one vector from either WH orbit 12 or 13 and the final 
vector from one of the orbits 23, 24, ...,28. As every normal vector in orbits 23 
to 28 are used, we find 6 x 36 = 216 triples. It is perhaps noteworthy that these 
latter six 'SIC-sensitive' WH orbits are precisely the Um invariant ones. 

In dimension 9, the SIC has 79,767 sets of nine linearly dependent vectors. 
Of these, 78,795 are predicted by Theorem 1 in Section 4. This large number of 
depedencies meant that we could not perform as complete computer searches as 
for the dimension 6 case. For example, we did not compute the scalar products 
between every pair of normal vectors, although we happen to know they form 
a basis (this is discussed in more detail in section 8). As with the dimension 6 
case, we can collect these 79,767 linearly dependent sets into orbits under the 
WH group. We find 987 WH orbits: 984 of length 81, 2 of length 27 and 1 of 
length 9. This gives 975 orbits invariant under the action of the Zauner unitary 
or one of its WH conjugates, with the remaining 12 invariant under Um (186 
orbits are invariant under both types). In odd dimensions where N = 3fc, M. 
takes the form 



When we repeat the calculations using a random vector from Hi, we again 
find exactly the same pattern of linearly dependent vectors. This suggests a dis- 
tinction between SICs in dimension 3 and SICs in higher dimensions divisible 
by 3. 'Special SICs' in dimension 3 gave rise to 12 sets of linear dependencies, 
while the others produced only 3. No SICs are 'special' in this sense in dimen- 
sions 6 and 9; instead we find that not only do all SICs have the same pattern 
of linearly dependent vectors, but this pattern can be be almost completely re- 
produced by acting with the WH group on any vector in the Zauner subspace 
(the pattern of orthogonalities among normal vectors cannot be reproduced). 
However, we did find one other instance of a SIC fiducial vector giving more 
linear dependencies than other vectors in the same subspace. In dimension 8, 
there is a SIC sitting in another eigenspace of the Zauner unitary, namely the 
one with the smallest dimension, W v 2. This SIC has more linear dependencies 
than a WH orbit starting with a random vector from W^: the SIC exhibits 
24,935,160 sets of eight linearly dependent vectors while a WH orbit using a 
vector in exhibits 24,756,984 (it is this latter value that is shown in Table 
2) [57]. SICs can always be found in this smaller subspace when the dimension 
equals 8 mod 9, which opens up the possibility of finding a connection between 
SICs and elliptic curves in this family of dimensions. 

We were unable to collect exhaustive information in dimension 12 regarding 
linear dependencies, so do not know how many sets of 12 linearly dependent 
vectors are in a SIC or a WH orbit starting from a random vector in Hi (or 
indeed any Zauner subspace). Similarly, we do not know whether their normal 




(33) 
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vectors exhibit any orthogonality structure analogous to the dimension 6 or 9 
cases. 

6 Embedded SICs in dimension N = 6 

In dimension N = 6, among 984 vectors normal to the linear dependent sets gen- 
erated from an arbitrary fiducial vector in the Zauner subspace (not necessarily 
a SIC fiducial), we found 30 sets of 4 normal vectors that form 2-dimensional 
SICs, i.e. in each set, the absolute values for the overlaps between the vectors are 
l/v3 and they lie on a 2-dimensional subspace. This phenomenon also happens 
in dimension N = 9, where 3-dimcnsional SICs are found among the normal 
vectors. Attempts to find such embedded SICs in dimension N = 12 yielded no 
positive result, perhaps only due to the impracticality of an exhaustive search. 
In this section, we give an explanation for the embedded SICs in dimension 
N = 6. In dimension N — 9, the phenomenon is not yet fully understood. 

The 2-dimensional embedded SIC sets found in dimension N = 6 are all 
of the form {|V) ,D 03 \tp) ,D 30 |V) ,D 33 \ip}}, where D^-'s are the displacement 
operators defined in Section 2. As will be proven in the following theorem, 
for such sets to form a 2-dimensional SIC, \tp) just needs to be any normalised 
eigenvector of Uz that is not in the Zauner subspace. Normal vectors of linearly 
dependent sets that contain three singlets in the case when the fiducial lies in 
the Zauner subspace just happen to satisfy this condition, as can been seen from 
Section 4. 

Theorem 4. In dimension N = 6, if is an eigenvector ofUz, then the 4 
vectors \ip) , Dq 3 | V) , D30 | V) and D 33 \tp) have constant absolute values for their 
pairwise overlaps. Furthermore, if |V) is not in the Zauner subspace, they span 
a 2-dimensional subspace and therefore constitute a SIC. 

Proof. If Uz IV) = A jV'), then \tp) = \U\ \ip)- The first part of the theorem can 
be proven by noticing that 

(VI D p |V) = (VI X*U z D p ulX |V> = (VI D Zp |V) (34) 

and that the Zauner symplectic matrix Z simply permutes the three points 03, 
30 and 33 in phase space according to the following table. 



p 


03 


30 


33 


Zp 


33 


03 


30 



To prove the second part of the theorem, we need to introduce the square 
root of the Zauner unitary. Let W be a symplectic matrix given by 

»-G -o 1 ) - 2 =(: =0-' <35 » 

Let £/yy be a unitary representative of W, with a phase chosen so that U^ = Uz- 
The eigenspaces of Uz and Uw are described in Table [4] where ui — e 27 ™/ 6 and 
r) = e 27 ™/ 3 = uj 2 . 
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Eigenspaces of Uz 


Hi 


'Zauner) 




T~L^2 


Eigenvalue 


1 


V 


2 


Dimensionality 


3 


2 


1 


Eigenspaces of Uy\> 


/Ci 


/C w 3 








Eigenvalue 


1 






CO 4 


w 2 


Dimensionality 


2 


1 


1 


1 


1 



Table 4: Eigenspace structure oi Uz and [Tyy. 



Next we define three new operators i?, S and T as follows. 
R={D 03 + D 30 + D 33 )/V3 

S = (£> 03 + ^ 2 A3o + w 4 £ 33 )/V3 (36) 
T = (Aw + oj 4 D 30 + oj 2 D 33 )/V3 

Note that ,D Q3 \i>) ,D 30 \iP) ,D 33 \i>)} and {\ip) ,R \ip) ,S\ip),T \ip)} have 
the same linear span. We will prove that \ip) e implies S\ip) = and 
R\tp) = \ip), and that G implies T \ip) = and R \tp) = - Therefore 
in both cases the span is 2-dimensional. 

The following properties of i?, S* and T can be derived straightforwardly from 
their definitions: 

S = T\ S 2 =T 2 = 0, R 2 = l, (37) 

ST = 1 + R, TS=1-R. (38) 

Moreover, ST/2 and TS/2 are rank-3 projection operators, and they are or- 
thogonal to each other. One can also verify the following commutation relations 
between R, S, T and C/yv 

U W R = RU W , U W S = lo 4 SU w , U W T = uj 2 TU w . (39) 

These equations tell us how R, S and T permute the eigenspaces of Uw- For 
example, if is an eigenvector of C/yv with eigenvalue u 2 , then UwS\ip) = 
ui 4 SUw \ip) — S\ip), so S \tp) is an eigenvector of C/yv with eigenvalue 1. The 
full action of S and T on the eigenspaces of [Tyy is described in the diagram 
below (R commutes with C/yy so it simply leaves the eigenspaces invariant). Let 
|fco), |fci), | fa), 1^3) and \ki) be non-zero eigenvectors of Uw belonging to the 
eigenspaces JCi, JC U , /C w 2, /C w 3 and /C w 4 respectively. Except for \ko), the rest 
of them are unique up to a scalar, because the corresponding eigenspaces are 
1-dimensional. As seen from the diagram, we have S \k±) = 0. We are going to 
prove that S |fc 4 ) = 0, so that S =0 for any e H v - 

Suppose otherwise, that S \k^) 7^ 0. It has to be the case that S\k ) = (for 
any choice of |fco)), because otherwise S |fco) will be a non-zero vector in /C w 4, 
which is 1-dimensional, which implies S |fco) = a ^4) for some non-zero a, and 
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n 




1~L\ 



>o 



Figure 2: The action of S (solid arrow) and T (dotted arrow) on the eigenspaces 
of U w - ICi (2D) and /C w 3 (ID) span Hi (3D). /C w 4 (ID) and K u (ID) span H v 
(2D). /C U 2 (ID) is H n 2. 

therefore S I&4) — a^S 2 \ko) — 0, contradicting the assumption that S \k4) ^ 0. 
Since S |/c4) is a non-zero vector in /C w 2, which is 1-dimensional, we must also 
have S \k4) = ft \k 2 ) for some non-zero (3, and therefore S \k 2 ) — fi~ 1 S 2 \k^) = 0. 
So from the assumption that S {k^} ^ 0, we have obtained = 5 |fci) = S |fc2) = 
S\k ), which implies = TS\ki) = TS\k 2 ) = TS\ko), which means that TS 
is orthogonal to the 4-dimensional subspace spanned by JCi, K,^ and /C w 2. This 
contradicts TS/2 being a rank-3 projection operator. 

Thus, we conclude that 5* \k^ = 0, and that S \ip) = for any \ip) € H v . The 
identity R\ip) = \ip) immediately follows from = TS\i/j) = (1 — R) \ip). Note 
that T \ip) (because TS is orthogonal to ST) is a non-zero vector in Hi. 
Therefore \ip), R\ip), S \ ip) and T \ indeed fully span a 2-dimensional subspace. 
It's a similar argument to show that T\k 2 ) = 0, which implies T\ip) — and 
R \ip) — — \ip) for any \ip) £ T-L^. This concludes the proof of the theorem. □ 

We can thus explain the occurrence of the 2-dimensional SICs one finds 
when N = 6. Unfortunately, one can't explain the 3-dimensional SICs one finds 
when TV — 9 in the same way. If one uses the construction described in the 
proof of Theorem 4 when N = 9 one obtains 9 vectors which do not belong 
to a 3-dimensional subspace and for which the overlaps are not constant (they 
actually take two different values). The construction also fails to produce SICs 
in dimensions 12 and 15. We summarize the situation in Table ?, which shows 
the dimensionality of the subspace spanned by the N 2 vectors obtained using 
the construction of Theorem 4 as a function of the eigenspace from which 
is chosen. 

Nevertheless, the fact remains that if, in dimension 9, instead of randomly 
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TV = 6 TV = 9 TV = 12 TV = 15 




4 
2 
2 



8 
7 
6 



12 

8 
8 



15 
15 
10 



Table 5: The dimension of the span of vectors created by the subgroup generated 
by D 3 and D 30 . 



choosing a vector in H v or H v 2 as in Theorem 4, one starts from a randomly 
chosen fiducial vector \ip) in the Zauner subspace (not necessarily a SIC fiducial), 
then one always does find 3-dimensional SICs among the normal vectors to the 
linearly dependent sets generated by | ■(/>). If this phenomenon repeated in higher 
dimensions it would open up the intriguing possibility, that one might be able 
to get a proof of SIC existence in this way. 

7 Generalisation to other symplectic unitaries 

We have been focussing so far on the Zauner unitary because, along with its 
conjugates, it is the symplectic unitary which is most relevant to the SIC exis- 
tence problem. However, there are many other symplectic unitaries which give 
rise to linear dependencies, as we now show. 

Before proceeding further we need to take account of a subtlety which is 
important in even dimensions. Suppose TV is even, so that TV — 2TV. Given 
G € SL(2,Zjv) let g r € SL(2,Zat) be Q reduced mod TV. We define 

1. The TV-order of Q to be its order considered as an element of SL(2, Z^-). 

2. The TV-order of Q to be the order of Q r considered as an element of 



This distinction was not important in the case of the Zauner matrix since its 
TV-order is the same as its TV-order. However, the two orders are, in general, 
different. Turning to Z^, we say that p, p' € Z^ are TV-distinct if p' ^ p mod 
TV. So as to be able to make statements which apply to both even and odd 
dimensions we will use the terms "TV-order" and "TV-distinct" in odd dimensions 
also, defining them in the obvious way. 

Let Q € SL(2, Zjy) (TV even or odd). We say that Q generates linear depen- 
dencies in all eigenspaces if, for each eigenvector \ip) of Ug, there exist pi, . . . p» 
which are TV-distinct and such that the vectors 



are linearly dependent. 

Continue to suppose Q e SL(2,Zjy), and let n be the TV-order of Q. We say 
that p is (/-full if the points p, Qp, . . . , G n ~ 1 p are TV-distinct. 

We then have 



SL(2,Z Ar ). 



D P1 M,...,D PW M 



(40) 
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Theorem 5. Let Q 6 SL(2,Z^) and let n be the N -order of Q . Suppose 
1. n > 1. 

<?. N is a multiple of n. 
3. Tr(U g ) ^0. 

^. There exist N points in 1? N which are N -distinct and Q-full 
Then Q generates linear dependencies in all eigenspaces. 

Proof. We first observe that Ug has at least one eigenspace with dimension 
< N/n. In fact, let 

n-1 

U g =e w J2° k Pk (41) 

fc=0 

be the spectral decomposition of Ug, where a = and e ie is a phase. By 
assumption 

n-1 

+ Tr([/ e ) = e ie ]T o~ k Tr(P fc ) (42) 

It follows that the numbers Tr(Pfc) cannot all be the same. Consequently we 
must have Tr(Pfc) < iV/n for at least one value of fe. 

It follows from condition[3]that it is possible to choose C*-full points pi , . . . , pjv/r, 
such that the TV points 

g r Pj , j = l,...,N/n, r = 0,...,n-l (43) 

are TV-distinct. Let \ip) be any eigenstate of Ug and define 

n-1 

l^fc) = $>- r,! A;. p >), j = l,...,N/n, fc = 0,...,n-l (44) 

r=0 

For all fc 

|Vi,fc>ef fc> j = l, — , N/n (45) 

where £ i , . . . , £ n are the eigenspaces t/p with the appropriate ordering (note 
that £k isn't necessarily the eigenspace corresponding to the projector P& defined 
above) . We showed that 

dim£ fc < N/n (46) 

for some k. For that value of k the vectors \ipi,k), ■ ■ ■ , iV'JV/n.fe) are linearly 
dependent, implying that the vectors 

Dg>~ Pj m, j = l,...,N/n, r = 0,...,n-l (47) 

are linearly dependent. □ 
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As an illustration of this theorem consider the parity matrix 



' ( o 1 iAS} 

We have n — 2 and Tr (U-p) = 2e t& (respectively e l8 ) if N is even (respectively 
odd). Also the number of £/-full points in 1? N is A 2 — 4 (respectively N 2 — 1) 
if N is even (respectively odd). So the theorem shows that V generates linear 
dependencies in every eigenspace for all even dimensions. 

Another application follows when considering the WH orbit MUBs men- 
tioned at the end of section 3. In every prime dimension N they are obtained 
by acting on a fiducial vector which is invariant under a Clifford unitary of order 
N. In this case our theorem can be used to explain why linear dependencies 
occur among the vectors in a MUB [26 . 

It will be seen that the conditions of the theorem are not very restrictive. 
There are A 2 points in 1? N so one would expect that the requirement that N of 
them be £/-full should not be difficult to satisfy. Similarly with the requirement 
that Tr([/g) 7^ 0. For a large class of dimensions we can make this statement 
sharp: 

Theorem 6. Let N — q\ . . .q u for u distinct odd prime numbers qi, . . . , q u . 
Then, with the sole exception of the identity matrix, every matrix e SL(2,Zjy) = 
SL(2,Zjv) whose order is a factor of N generates linear dependencies in every 
eigenspace. 

The proof of this theorem will occupy us for the rest of this section. It 
depends on the tensor product representation of the Clifford group which is 
described in Appendix B in |13j . and which we now briefly review. For each 
x G Zjv and j in the range 1 , . . . , u let Xj be the unique element of 1 q . such 
that Xj = x mod qj. Then the Chinese Remainder Theorem |2S] tells us that 
the map 

x (xi,...,x u ) (49) 
is an isomorphism of the ring Zjv onto the ring Z qi x • • ■ x Z Qu with inverse 

NkiXi Nk u x u 
(xi,...,x u ) i-> 1 1 (50) 

where Kj is the multiplicative inverse of N/qj considered as an element of "L qj . 
Given 

Pl ^ez% g=( a f)eSL(2,Zjv) (51) 



Pi) V7 5 



define 
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Then define 

<j>: Z% -S-Z^ 0(p) = Pl ® ■ ■ ■ ® p„ (53) 

V>: SL(2,Zat) -> SL(2,Z gi )<»••• ®SL(2,Z,J ^(0) = Qx ® • • • (8 0« (54) 

The map <f> (respectively t/;) is an isomorphism of modules (respectively groups). 
We have 

<P(Qp) = ^(0)«/>(p) (55) 
for all G, p. Let T-L n be n-dimensional Hilbert space and let 

V: Hn ~> H qi (&■■■<£) Hg u (56) 

be the unitary defined by 

V\x) N = \xi) qi ®---®\x u ) q „ (57) 

where |0)„, . . . , \n — 1)„ is the standard basis in dimension n. Then 

VD p V^ = D P1 <8> • • • ® D p „ (58) 
FZ7 e Ft = Z7 0l ® • • • <g) U Qu (59) 

modulo the phase ambiguity in the definitions of C/g, Ug 1 , . . . , i/c; u . 

Lemma 7. Let N = q± . . . q u for u distinct odd prime numbers qi, . . . , q u . Let 
Q G SL(2,Zjv) &e arbitrary. Then the number of Q- full points in 1j 2 n is > 
JV(gi-l)...( ff „-l). 

Proof. Let if be the number of C?-full points in 1? N and -fTj the number of Qj- 
full points in 1? . Use the isomorphisms defined in Equations ([53")) and ([54")) to 
identify with Q\® ■ ■ ■ ® Q u and p with Pl <£> • • • <g> p u . Let n be the order of Q 
and rij the order of . Then 

n = LCM(m,...,n u ) (60) 

Suppose pj is C/j-full in Z^. for each j. Then 

g r p = p (61) 

implies r is a multiple of nj for each j, which in turn implies r is a multiple of n. 
So p is (/-full in 1? . So K > K\ . . . K u . The problem thus reduces to showing 
that Kj > qj(qj — 1) for all j. 

To see this we use the analysis in Let tj — Tr(Qj) and let Qj (respec- 
tively Nj) be the set of quadratic residues (respectively non-residues) in Z qj . If 
tj-4g Qj then Qj is conjugate to a diagonal matrix 



a 
a" 1 



(62) 
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If pj is non-zero then 



SjPi = Pj 



(63) 



if and only if r is a multiple of nj. So Kj = qj — 1. If i| — 4 S Nj we 
can diagonalize CJj over the extension field F ? 2 . So we again have that Qj is 
conjugate to a matrix of the form 



a 
a' 1 



(64) 



where a now belongs to F ? 2 instead of Z 9:j . We again deduce that Kj = qj — 1. 
Finally, if = ±2 we have that Qj is conjugate to the matrix 



±1 /3 
±1 



for some (3. If /3 = then Kj = qj — 1; otherwise ifj — qj — qj. 



(65) 



□ 



Lemma 8. Let N = qi . . . q u for u distinct odd prime numbers qi, . . . , q u . Then 
\Tr(Ug)\ > 1 for all Q € SL(2,Zat). 

Proof. In view of Equation (fST))) it is enough to show that 

|Tr([/ Sj )| >1 



for all j. To see this let 



a (3 
7 S 



We have 129 



3 1 „ifl _o7a: 2 



Z^£C = 



ax) (a; | 



0^0 
/3 = 



(66) 
(67) 

(68) 



where t, = e 9 j and is a phase. So 



Tr(U gi ) = {e ie 



2- 1 /3- 1 (t-2)a 
*=0 W j 



/3^0 

/3 = 0, a ^ 1 
/3 = 0, a = 1 



(69) 



where t — Tr(C^ ) and Uj — e ^ . Performing the Gauss sums |28j we find 



Tr(C/e 3 ) = 



x]e i6 ^fq-j 
e w 



/3 ^ 0, i = 2 

P = 0, a^l 

= 0, a = 1, 7^0 

= 0, a= 1, 7 = 



(70) 
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where 




1 qj = 1 mod 4 

i q 3 ; = — 1 mod 4 



where (— ) is the Legendre symbol, and where we have used the multiplicative 

property of the Legendre symbol together with the fact that ) = (|-) for all 
non-zero x. We see that in every case 



1^)1 = 1. VQJ ot 1j (72) 



□ 



Proof of Theorem [3J Immediate consequence of Theorem [5] and Lemmas [7] |5] 



It would be interesting to see if one can generalise Theorem [6] to the case of 
even dimensions, and odd dimensions for which the prime factorization contains 
primes occurring with multiplicity > 1. 

8 Linear dependencies in the monomial represen- 
tation 

In the motivating three-dimensional example there was an interesting interplay 
between linear dependencies in WH orbits on the one hand, and preferred or- 
thonormal bases on the other. It is not difficult to find fiducial vectors giving 
rise to a WH orbit with N linearly dependent sets of N vectors each, in such a 
way that their normal vectors form an orthonormal basis: if the generator Z is 
diagonal, any fiducial vector with a zero entry will give rise to such a linearly 
dependent set under repeated action with Z . Acting with X on the fiducial will 
permute its entries, and one ends up with N linearly dependent sets in the orbit. 
However, these fiducials are not very interesting, and these orbits certainly do 
not include SICs. 

The situation is more interesting if the dimension is a square number, N = 
n 2 . For this case it has recently been shown that there exists a choice of basis 
in which the entire Clifford group is represented by monomial matrices, i.e. by 
permutation matrices in which the unit entries have been replaced by phase 
factors |13| . This representation stems from the observation that in square 
dimensions 



Thus X n and Z n can be simultaneously diagonalised, and the resulting basis 
turns out to have the above property. Denote the basis vectors by \r, s), where 
r, s are integers modulo n, and let us refer to the basis as the monomial basis. 
The monomial representation of the Clifford group is given by 



□ 



[X n ,Z n ] = 



(73) 




q r \r, 0) if s + 1 = mod n 



|r, s + 1) if s + 1 7^ mod n 



(74) 
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Z\r,s) = w s |r- l,s) 



(75) 



U G \r,s) =e i0 T f} W 2 -2 S s'W)| 5r _ 7S + mi s,-p r + as + map) (76) 
where 

ia f n is odd 

s = —pr + as + ma m={ „ . (77) 

[ n is even v y 

Note that Equation ([75)) assumes /3 is relatively prime to TV, otherwise one has 
to use the decomposition given in Equation (fT4]) of Section 2. Evidently the N 
orthogonal rays that make up the basis is in itself an orbit of the full Clifford 
group. Moreover, one sees by inspection that the vector |0, 0) is left invariant 
by the full symplectic group whenever N is odd. In this sense the monomial 
basis is a distant cousin of the 'special SIC in three dimensions — or perhaps of 
the MUB in three dimensions. 

Now suppose that 3|n, and consider the Zauner matrix 

*-(!:!)■ ™ 

If we scrutinise the explicit expression (I76|) we find that there are three diag- 
onal elements in the unitary operator Uz- If we adjust the phase e according 
to Zauner's recipe one of them equals rj, and the other two equal 1. It follows 
that any vector in the Hi eigenspace must have one component equal to 0. Since 
the WH group contains N diagonal elements (namely the elements of the sub- 
group generated by X n and Z n ) this means that the WH orbit generated from 
a fiducial vector in Hi contains N vectors confined to a subspace orthogonal to 
one of the vectors in the monomial basis. Since the entire WH group is given 
by monomial matrices the orbit can be divided into N such linearly dependent 
sets, spanning N subspaces each orthogonal to a vector in the monomial basis. 
In particular this applies to the known SICs in these dimensions. 

A similar argument applies if 2\n. We consider the symplectic matrix 

*-(-?;) <™> 

After adjusting the arbitrary phase so that Up = 1 and so that the largest 
eigenspace corresponds to the eigenvalue 1, we find that Ujr has two diagonal 
non-zero elements, one equal to 1 and the other equal to uj n ^ 4 . Following 
the same logic as above we see that every WH orbit that includes a fiducial 
invariant under Ujr must be divided into N linearly dependent sets of TV vectors 
orthogonal to one of the vectors in the monomial basis. 

Variants of this argument can be used also if the dimension N = kn 2 , where 
k is any integer. In this case Hilbert space splits into a direct sum of n 2 copies 
of a fc-dimensional Hilbert space, and there exists a representation of the Clif- 
ford group in which every group element acts within and permutes these sub- 
spaces [30 . This representation is fc-nomial, in the sense that it is given by a 
permutation matrix of size n 2 in which the unit entries have been replaced by 
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matrices of size k. For instance, one finds in dimension eight that any vector in 
the smallest eigenspace of the Zauner unitary has two zero entries, and again 
there will be linear dependencies in a WH orbit arising from a fiducial in this 
eigenspace. A SIC fiducial does in fact exist in this eigenspace |12| . 

9 Conclusion 

The Hesse configuration singles out 'special SICs' in dimension 3 by the higher 
number of linear dependencies among their vectors. We searched for similar 
patterns in higher dimensions through exhaustive numerical searches for sets of 
N linearly dependendent vectors in SICs in dimensions TV = 4 to 8 and partial 
searches in dimensions N — 9 and 12. In dimension 8, the SICs with fiducials in 
the smallest subspace of the Zauner unitary, H v 2 , also stand out for exhibiting 
a higher number of linear dependencies than other WH orbits with a fiducial 
vector from the same subspace. This could extend the dimension 3 connection 
between elliptic curves and SICs (which gives rise to the Hesse configuration) 
to dimension 8 and, as SIC fiducials always lie in H v 2 for dimensions equal to 8 
mod 9, possibly to an infinte family of dimensions. 

In the dimensions we looked at, sets of N SIC vectors are usually linearly in- 
dependent. We proved that it is always possible to find dependencies in SICs for 
dimensions divisible by 3 (where the SIC fiducials sit in "Hi) and for dimensions 
equal to 8 mod 9 (where the SIC fiducials sit in H^). However, this dependency 
structure in dimensions divisible by 3 also occurs for any WH orbit when the 
fiducial vector is taken from Hi and so is not a property uniquely tied to SICs. 
We also proved that linear dependencies will arise for WH orbits whose fiducial 
vector is in Hi for dimensions N = 3k; WH orbits whose fiducial vector is in H v 
for dimensions N — 3k and N = 3k + 1; and WH orbits whose fiducial vector 
is in H v 2 for all dimensions. These theorems do not account for all the linear 
dependencies; our numerical results often revealed a higher number of linearly 
dependent sets of vectors that the theorems predict. 

Looking in detail at linear dependencies in dimension 6 revealed some sur- 
prising structure. Firstly, the normal vectors (i.e. the vectors orthogonal to each 
5-dimcnsional subspace spanned by 6 linearly dependent vectors) collect into or- 
thogonal sets of 4 when the HW orbit is a SIC. When the orbit is not a SIC, these 
normal vectors lose their orthogonality pattern. Secondly, the normal vectors of 
dependencies from both SICs and non-SICs form 30 2-dimensional SICs. This 
repeats in dimension 9, where some of the normal vectors form 3-dimensional 
SICs. We did not find 4-dimensional SICs from normal vectors in dimension 12, 
but our search was not exhaustive and we cannot rule out their existence either. 

We also proved a number of theorems for dependencies using unitaries other 
than the Zauner unitary. Our results represent the beginning of a theory for 
when linear dependencies occur in WH orbits. This is of interest in itself. Also, 
just possibly, they may contain an important clue for the SIC existence problem. 
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